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depth- sensing  indentation  or  nanoindenta¬ 
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the  mechanical  response  of  materials  from 
metals  and  ceramics  to  polymeric  and 
biological  materials.  The  additional  levels 
of  control,  sensitivity,  and  data  acquisition 
offered  by  instrumented  indentation 
systems  have  resulted  in  numerous 
advances  in  materials  science,  particularly 
regarding  fundamental  mechanisms  of 
mechanical  behavior  at  micrometer  and 
even  sub -micrometer  length  scales. 
Continued  improvements  of  instrumented 
indentation  testing  towards  absolute 
quantification  of  a  wide  range  of  material 
properties  and  behavior  will  require 
advances  in  instrument  calibration, 
measurement  protocols,  and  analysis  tools 
and  techniques.  In  this  paper,  an  overview 


of  instrumented  indentation  is  given  with 
regard  to  current  instrument  technology 
and  analysis  methods.  Research  efforts  at 
the  National  Institute  of  Standards  and 
Technology  (NIST)  aimed  at  improving 
the  related  measurement  science  are 
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1.  Introduction 

Indentation  or  hardness  testing  has  long  been  used 
for  characterization  and  quality  control  of  materials, 
but  the  results  are  not  absolute  and  depend  on  the  test 
method.  In  general,  traditional  hardness  tests  consist  of 
the  application  of  a  single  static  force  and  correspon¬ 
ding  dwell  time  with  a  specified  tip  shape  and  tip 
material,  resulting  in  a  hardness  impression  that  has 
dimensions  on  the  order  of  millimeters.  The  output  of 
these  hardness  testers  is  typically  a  single  indentation 
hardness  value  that  is  a  measure  of  the  relative  pene¬ 
tration  depth  of  the  indentation  tip  into  the  sample.  For 


1  Current  affiliation:  Army  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD  21005. 


example,  the  Rockwell  hardness  scales  are  distin¬ 
guished  by  the  amount  of  static  force  applied  (e.g., 
100  kg  for  Rockwell  B  compared  to  150  kg  for 
Rockwell  C);  the  tip  shape  (e.g.,  a  spherical  tip  with  a 
diameter  of  1.6  mm  for  Rockwell  B  compared  to  a 
conical  tip  that  has  a  spherical  apex  with  a  radius  of 
200  pm  for  Rockwell  C);  and  the  tip  material  (e.g.,  steel 
for  Rockwell  B  compared  to  diamond  for  Rockwell  C) 
[1].  Harder  metals  are  more  appropriately  characterized 
using  higher  forces  and/or  pyramidal  tips,  whereas 
lower  forces  and  spherical  tips  are  used  on  softer 
metals.  Similarly,  different  durometers,  which  are  used 
to  characterize  the  mechanical  resistance  of  polymers, 
have  different  spring  constants  and  either  a  flat  conical 
tip,  a  sharp  conical  tip,  or  a  spherical  tip,  as  specified 
in  ASTM  D  2240  Standard  Test  Method  for  Rubber 
Property — Durometer  Hardness.  Durometers  with 
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higher  spring  constants  and  sharp  tips  can  be  used  to 
evaluate  stiffer  polymers  compared  to  durometers  with 
lower  spring  constants  and/or  flat  or  spherical  tips. 

Instrumented  indentation,  also  known  as  depth¬ 
sensing  indentation  or  nanoindentation,  is  increasingly 
being  used  to  probe  the  mechanical  response  of  materi¬ 
als  from  metals  and  ceramics  to  polymeric  and  bio¬ 
logical  materials.  In  contrast  to  traditional  hardness 
testers,  instrumented  indentation  systems  allow  the 
application  of  a  specified  force  or  displacement  history, 
such  that  force,  P,  and  the  displacement,  h,  are  con¬ 
trolled  and/or  measured  simultaneously  and  continu¬ 
ously  over  a  complete  loading  cycle.  Additionally,  the 
extremely  small  force  and  displacement  resolutions, 
often  as  low  as  -1  pN  and  ~0.2  nm,  respectively,  or 
lower  for  some  systems,  are  combined  with  very  large 
ranges  of  applied  forces  and  displacements  (tens  of  pN 
to  hundreds  of  mN  or  larger  in  force  and  tens  of  nm  to 
tens  of  pm  or  larger  in  displacement)  to  allow  a  single 
instrument  to  be  used  to  characterize  nearly  all  types  of 
material  systems.  Recent  developments  include 
improved  user  control  over  loading  histories,  improved 
system  and  software  robustness  allowing  much  higher 
levels  of  test  automation,  and  the  use  of  dynamic 
oscillation  for  improved  sensitivity  and  additional 
testing  capabilities. 

The  additional  levels  of  control,  sensitivity,  and 
data  acquisition  offered  by  instrumented  indentation 
systems  have  resulted  in  numerous  advances  in  materi¬ 
als  science,  particularly  regarding  fundamental  mecha¬ 
nisms  of  mechanical  behavior  at  micrometer  and  even 
sub-micrometer  length  scales.  Instrumented  indentation 
systems  have  been  used  to  study,  for  example,  disloca¬ 
tion  behavior  in  metals  [2-4],  fracture  behavior  in 
ceramics  [5-6],  mechanical  behavior  of  thin  films 
[7-10]  and  bone  [11],  residual  stresses  [12],  and  time- 
dependent  behavior  in  soft  metals  [13-15]  and  poly¬ 
mers  [16-22].  Additionally,  lateral  probe  motion  is  also 
being  incorporated  to  explore  tribological  behavior  of 
surfaces,  including  scratch  resistance  of  coatings 
[23,  24]  and  wear  resistance  of  metals  [25].  The  extent 
to  which  these  techniques  can  be  used  to  quantify 
macroscopic  material  properties  is  continually  being 
explored.  Indentation  values  of  elastic  modulus  are 
routinely  calculated  and  compared  to  values  from 
macroscopic  mechanical  tests.  Further,  a  number  of 
researchers  are  utilizing  techniques  such  as  dimension¬ 
al  analysis  and  finite  element  modeling  to  relate  mate¬ 
rial  constitutive  behavior  to  indentation  force-displace¬ 
ment  data,  and  even  to  predict  stress- strain  behavior 
based  on  instrumented  indentation  data  [26,  27]. 


Continued  improvements  of  instrumented  indenta¬ 
tion  testing  towards  absolute  quantification  and 
standardization  for  a  wide  range  of  material  properties 
and  behavior  will  require  advances  in  instrument 
calibration,  measurement  protocols,  and  analysis  tools 
and  techniques.  In  this  paper,  an  overview  of  instru¬ 
mented  indentation  is  given  with  regard  to  current 
instrument  technology  and  analysis  methods.  Research 
efforts  at  the  National  Institute  of  Standards  and 
Technology  (NIST)  are  then  discussed,  which  are 
aimed  at  improving  the  related  measurement  science. 

2.  Overview  of  Instrumented  Indentation 

2.1  Instrumentation 

Many  instrumented  indentation  systems  can  be 
generalized  in  terms  of  the  schematic  illustration  shown 
in  Fig.  1  [28,29].  Commercially  available  systems 
include  those  developed  by  Nano  Instruments1  (Oak 
Ridge,  TN;  now  part  of  MTS  Systems  Corp.)  [28], 
Hysitron,  Inc.  (Minneapolis,  MN)  [30],  Micro 
Materials  Ltd.  (UK)  [31],  CSIRO  (Australia)  [32],  and 
CSEM  (now  CSM)  Instruments  (Switzerland)  [23]. 
Force  is  often  applied  using  either  electromagnetic  or 
electrostatic  actuation,  and  a  capacitive  sensor  is  typi¬ 
cally  used  to  measure  displacement.  In  the  MTS  Nano 
Instruments,  Hysitron,  CSIRO,  and  CSM  Instruments 
systems,  the  axis  of  the  indenter  is  vertical  (see  Fig.  1), 
while  in  the  Micro  Materials  system,  it  is  horizontal. 
Also,  the  MTS  Nano  Instruments,  Micro  Materials, 
CSIRO,  and  CSM  Instruments  systems  apply  force  and 
measure  displacement  through  separate  means,  while 
the  Hysitron  system  uses  the  same  transducer  for  both 
force  application  and  displacement  measurement. 
Regardless  of  the  operational  differences,  the  raw  force 
and  raw  displacement  data  for  these  systems  are  always 
coupled  due  to  the  leaf  springs,  and  in  fact,  a  plot  of  the 
raw  force  as  a  function  of  the  raw  displacement  with 
the  tip  out  of  contact  is  characteristic  of  the  spring 
stiffness.  In  the  following  sections,  the  measurement 
methods,  analysis  techniques,  and  calibration  issues 
presented  are  applicable,  in  general,  to  current  com¬ 
mercial  instrumented  indentation  systems. 

1  Certain  commercial  equipment,  instruments,  or  materials  are  iden¬ 
tified  in  this  paper  to  foster  understanding.  Such  identification  does 
not  imply  recommendation  or  endorsement  by  the  National  Institute 
of  Standards  and  Technology,  nor  does  it  imply  that  the  materials 
or  equipment  identified  are  necessarily  the  best  available  for  the 
purpose. 
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Fig.  1.  Schematic  illustration  of  an  instrumented  indentation  system. 


2.2  Measurement  Methods 


Some  of  the  earliest  research  related  to  the  develop¬ 
ment  of  instrumented  indentation  is  found  in  the  early  - 
to  mid-1980's  [29,  31,  33,  34].  Since  that  time,  the  con¬ 
tinuing  evolution  of  computer  technology  has  led  to 
improved  control  over  instrumented  indentation  testing. 
Although  some  specialized  research  instruments  are 
displacement  controlled  (for  example,  the  inter-facial 
force  microscope  [35]),  most  indentation  systems  are 
force-driven  devices,  such  that  the  force  history  is  most 
easily  controlled  with  displacement  control  achieved 
using  signal  feedback.  Typical  force-time 
( P-t )  and  force-displacement  (P-h)  data  are  shown  in 
Fig.  2.  In  Fig.  2a  and  Fig.  2b,  indentation  data  taken 
using  a  constant  loading  rate,  P  ,  is  shown  for  poly 
(methyl  methacrylate),  whereas  in  Fig.  2c  and  Fig.  2d, 
loading  data  are  shown  for  constant  P/P  for  poly¬ 
styrene.  In  both  cases,  the  loading  segment  is  followed 
by  a  hold  segment,  which  is  then  followed  by  an 
unloading  segment  using  a  constant  unloading  rate.  In 
Fig.  2c  and  Fig.  2d,  a  hold  period  near  the  end  of  the 
unloading  segment  is  also  shown.  Assuming  the 
mechanical  behavior  of  the  material  is  not  a  function 
of  time,  the  ratio  of  loading  rate,  P  ,  to  force,  P,  has 
been  related  to  a  proposed  expression  for  strain  rate  for 
an  indentation  measurement,  £  ,  through  the  follow¬ 
ing  relation  [36]: 


h~  p  P  H 
^  / 


(i) 


Here,  H  =  P/A  is  the  hardness,  H  =  d H  /  d t  is  the  time 
rate  of  change  during  testing  of  the  hardness,  and  /3 
describes  the  shape  of  an  idealized  indentation  tip, 
where  the  cross-section  area,  A,  is  related  to  the 


distance,  h,  from  the  tip  apex  by  A  =  ch^.  Note  that  in 
Eq.  (1),  /3  cannot  be  0,  so  that  the  right  side  of  this 
equation  is  not  valid  for  flat  punch  geometry.  For  a  tip- 
sample  combination  for  which  H  is  constant  with  depth, 
constant  P/P  testing  yields  a  constant  indentation 
strain  rate.  For  many  materials,  H  is  a  function  of  £  , 
which  can  then  be  characterized  using  these  types  of 
tests.  Further,  because  H  =  P/A  represents  a  mean  pres¬ 
sure  or  stress  during  indentation,  constant  PIP  testing 
could  be  used,  for  example,  to  study  creep  behavior 
under  “constant  H”  (or  equivalently,  “constant  mean 
pressure”  or  “constant  mean  stress”)  conditions  in  time- 
dependent  materials. 

A  more  ubiquitous  test  method  used  in  force- 
controlled  instrumented  indentation  employs  the  use  of 
a  constant  loading  rate,  or  for  a  displacement-controlled 
system,  a  constant  displacement  rate.  In  either  case,  no 
feedback  is  necessary,  which  lends  to  the  simplicity  of 
these  methods.  Of  course,  feedback  could  be  used  to 
produce  a  constant  displacement  rate  in  a  force- 
controlled  system  or  a  constant  loading  rate  in  a 
displacement-controlled  system.  The  former  approach 
is  perhaps  more  useful,  because  displacement,  h, 
and  displacement  rate,  h ,  are  linked  to  indentation 
strain  and  strain  rate,  whereas  force  and  loading  rate  are 
more  difficult  to  link  to  either  stress  or  strain.  For  exam¬ 
ple,  for  conical  or  pyramidal  tip  geometry,  a  nominal 
indentation  strain  is  related  to  the  characteristic  includ¬ 
ed  angle  or  angles  of  the  tip.  For  a  paraboloidal  tip, 
indentation  strain  is  related  to  the  ratio  of  the  contact 
radius  to  the  tip  radius,  and  the  contact  radius  is  direct¬ 
ly  related  to  h  [16,  21].  For  any  tip  geometry,  the  inden¬ 
tation  strain  rate  can  be  calculated  from  the  ratio  h/h, 
as  in  Eq.  (1).  However,  the  mean  stress  or  hardness,  H , 
is  the  ratio  of  force,  P,  to  contact  area,  A,  which  in  gen¬ 
eral  is  related  to  displacement  through  the  tip  geometry. 
Only  in  the  case  of  a  flat  punch,  where  A  is  constant 
with  h,  is  H  a  function  of  force  only.  As  shown  in 
Eq.  (1),  P/P  can  be  related  to  indentation  strain  rate, 
but  only  through  additional  calculations  of  H  and  H . 

Other  quasi- static  test  methods  that  have  been  used  to 
characterize  time-dependent  response  are  indentation 
creep  tests  [13-15,  17,  21]  and  indentation  stress  relax¬ 
ation  tests  [21].  As  in  a  traditional  tensile  creep  test, 
force  is  held  constant  in  an  indentation  creep  test  and 
displacement  or  strain  is  monitored  [1].  However,  for 
the  case  of  tensile  creep,  the  change  in  cross-section 
area  as  the  sample  elongates  is  typically  small  during 
most  of  the  test,  such  that  constant  force  is  essentially 
equivalent  to  constant  stress.  In  the  indentation  creep 
test,  the  gradual  displacement  of  the  tip  into  the  sample 
causes  a  substantial  change  in  the  contact  area,  as 
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Fig.  2.  Force-time  (a  and  c)  and  corresponding  force-displacement  (b  and  d)  indentation  data:  (a  and  b)  data  taken  on  poly(methyl  methacrylate) 
using  a  constant  P  loading;  (c  and  d)  data  taken  on  polystyrene  using  constant  PIP  loading.  Both  data  sets  include  a  10  s  hold  period  between 
loading  and  unloading,  and  the  data  in  (c)  and  (d)  include  an  additional  hold  period  near  the  end  of  the  unloading  segment. 


shown  in  Fig.  3,  such  that  both  the  stress  and  the  strain 
fields  evolve  during  the  test,  complicating  analysis  and 
comparison  to  bulk  measurements.  For  the  indentation 
stress  relaxation  test,  displacement  is  held  constant 
(through  feedback  in  a  force-controlled  system)  and  the 
gradual  decrease  in  force  is  recorded.  Because  dis¬ 
placement  is  closely  linked  to  strain  in  an  indentation 
experiment,  constant  indentation  displacement  gives 
constant  indentation  strain.  Thus,  data  from  this  type  of 
measurement  method  are  more  easily  analyzed  and 
compared  to  traditional  bulk  stress  relaxation  data. 


However,  because  very  fast  initial  force  and  displace 
ment  rates  are  normally  used  in  creep  and  stress  relax¬ 
ation  tests,  respectively,  the  use  of  feedback  combined 
with  instrument  limitations  typically  causes  issues 
regarding  how  fast  the  force  or  displacement  can  be 
applied  and  then  held  constant  without  substantial  over¬ 
shooting  or  undershooting,  which  can  affect  the  quality 
particularly  of  the  short-time  data.  An  example  of 
indentation  stress  relaxation  data  taken  on  a  force- 
controlled  system  is  shown  in  Fig.  4  to  illustrate  this 
problem. 
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Fig.  3.  Force-time  (a)  and  corresponding  displacement-time  (b)  and  contact  area-time  data  (c)  for 
an  indentation  creep  test  using  a  force-controlled  system. 
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Fig.  4.  Displacement-time  (a)  and  corresponding  force-time  (b)  for  an  indentation  stress-relaxation 
test  using  a  force-controlled  system  and  feedback  to  control  displacement. 


Fig.  5.  Schematic  illustration  of  a  dynamic  model  for  an  instru¬ 
mented  indentation  system.  Kf  represents  the  load-frame  stiffness, 
Ks  represents  the  stiffness  of  the  springs,  D  and  mi  represent 
the  damping  characteristics  and  mass,  respectively,  of  the  instru¬ 
ment,  and  S  and  C  represent  the  storage  and  loss  components,  re¬ 
spectively,  of  the  mechanical  impedance  related  to  the  tip- sample 
contact. 


Recently,  test  methods  have  been  used  that  combine 
dynamic  oscillation  with  the  quasi-static  testing  capa¬ 
bilities  of  instrumented  indentation  systems.  A  typical 
dynamic  model  of  the  indentation  system  represented 
in  Fig.  1  is  shown  schematically  in  Fig.  5.  When 
dynamic  oscillation  is  applied,  it  is  most  often  super¬ 
posed  over  a  quasi-static  force  history  using  a  small 
force  or  displacement  amplitude  and  a  frequency  in  the 
range  of  1  Hz  to  300  Hz  [19,  37].  As  will  be  discussed 
further  in  Sec.  2.3,  the  equations  developed  for  the 
dynamic  model  are  used  to  determine  the  contact  stiff¬ 
ness  throughout  the  force  history,  which  in  turn  allows 
the  contact  area  and  sample  modulus  to  be  estimated 
throughout  the  force  history.  This  technique  is  thus 
particularly  useful  for  characterizing  modulus  as  a 
function  of  depth,  estimating  changes  in  contact  area 
during  an  indentation  creep  test  (see  Fig.  3c),  and 
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exploring  the  storage  and  loss  responses  of  materials, 
for  example.  Additionally,  the  statistical  sampling  of 
indentation  testing  is  dramatically  improved  over  get¬ 
ting  one  set  of  values  for  a  given  test  from  analysis  of 
the  unloading  curve  slope  (see  Sec.  2.3).  Also,  certain 
parameters  related  to  the  system  dynamics  will  change 
dramatically  at  the  onset  of  tip- sample  contact,  signifi¬ 
cantly  improving  the  ability  to  determine  this  point  in 
the  indentation  data. 


2.3  Analysis  Techniques 

The  analysis  of  force-displacement  curves  produced 
by  instrumented  indentation  systems  is  often  based  on 
work  by  Doerner  and  Nix  [29]  and  Oliver  and  Pharr 
[28].  Their  analyses  were  in  turn  based  upon  relation¬ 
ships  developed  by  Sneddon  [38]  for  the  penetration  of 
a  flat  elastic  half  space  by  different  probes  with  partic¬ 
ular  axisymmetric  shapes  (e.g.,  a  flat-ended  cylindrical 
punch,  a  paraboloid  of  revolution,  and  a  cone).  These 
elasticity-based  analyses  are  normally  applied  to  the 
unloading  data  of  an  indentation  measurement,  assum¬ 
ing  the  unloading  behavior  of  the  material  is  character¬ 
ized  by  elastic  recovery  only.  In  general,  the  relation¬ 
ships  between  penetration  depth,  h ,  and  force,  P,  during 
unloading  can  be  represented  in  the  form 

P  =  a(h-hf)m  (2) 

The  parameter  a  contains  geometric  constants,  the 
sample  elastic  modulus,  E ,  the  sample  Poisson's  ratio, 
v,  the  indenter  elastic  modulus,  E{,  and  the  indenter 
Poisson's  ratio,  v-;  the  parameter  hf  is  the  final  unload¬ 
ing  depth;  and  m  is  a  power  law  exponent  that  is  relat¬ 
ed  to  the  geometry  of  the  indenter  (see  Table  1).  A  non¬ 
linear  power  law  fit  to  the  unloading  data,  where  a,  hf , 
and  m  are  fitting  parameters,  often  yields  a  good 
estimate  of  the  data,  as  does  a  smooth  spline  fit  [20]. 
However,  previous  research  at  NIST  has  shown  that  the 
practice  of  arbitrarily  fitting  portions  of  the  unloading 
data  can  introduce  bias  into  the  calculations,  and  that 
data  should  be  removed  for  curve  fitting  purposes  only 
when  the  corresponding  residual  errors  do  not  conform 
to  the  assumptions  underlying  the  curve  fitting  methods 
[20].  Once  an  appropriate  fit  is  obtained,  a  derivative, 
dP/dh ,  applied  at  the  maximum  loading  point  (/zmax, 
Pmax)  should  yield  information  about  the  state  of  contact 
at  that  point.  This  derivative  is  termed  the  contact  stiff¬ 
ness,  S ,  and  is  given  analytically  by 


S  =  2aEr=^=EI4A. 

V7T 


Table  1.  Theoretical  values  of  parameters  m  and  £,  both  of  which  are 
related  to  the  contact  geometry,  for  three  axisymmetric  tip  shapes 
[38],  where  m  is  the  power  law  exponent  of  Eq.  (2)  and  £  is  a  factor 
used  in  determining  the  contact  depth  [see  Eq.  (6)]. 


Tip  geometry 

m 

£ 

Flat-ended  cylindrical  punch 

1 

1 

Paraboloid  of  revolution 

1.5 

0.75 

Cone 

2 

2  (tu— 2)/tc 

In  this  equation,  the  cross  section  of  the  indenter  is 
assumed  to  be  circular  in  relating  the  contact  radius,  a , 
to  the  projected  area  of  tip- sample  contact,  A.  A  small 
correction  is  sometimes  applied  for  non-circular  cross 
sections  [26]  and  other  corrections  have  been  suggest¬ 
ed  [39].  The  reduced  modulus,  ET,  accounts  for  elastic 
deformation  of  both  the  indenter  and  the  sample  and  is 
given  by 


1  _  (1  —  y2)  |  (1-v*) 

Ex  E  Ex 


(4) 


In  Fig.  6,  an  indentation  force-displacement  curve  is 
illustrated  along  with  several  important  parameters 
used  in  the  Oliver  and  Pharr  analysis.  The  measured 
stiffness,  S *,  is  the  slope  of  the  tangent  line  to  the 
unloading  curve  at  the  maximum  loading  point  (Pmax) 
and  is  given  by 

=  am(Ks  ~  K  T'  ■  (5) 


Displacement,  h 


hmax  =  maximum  displacement 
hf  =  final  depth 


(3) 
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Fig.  6.  An  indentation  force-displacement  curve  in  which  several 
important  parameters  used  in  the  Oliver  and  Pharr  analysis  are 
illustrated. 
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The  parenthetic  subscript  denotes  that  the  derivative  is 
evaluated  at  the  maximum  loading  point.  Physically, 
the  elastic  recovery  of  the  material  is  instantaneous 
upon  unloading,  so  the  true  elastic  response  of  the 
material  can  only  be  evaluated  at  t-  0+,  where  unload¬ 
ing  initiates  at  t  -  0.  When  the  displacement,  h,  is  the 
total  measured  displacement  of  the  system,  S  *  is  the 
total  system  stiffness.  After  successful  calibration  of 
the  load-frame  compliance  (see  Sec.  2.4),  the  displace¬ 
ment  of  the  load  frame  is  removed  so  that  h  represents 
only  the  displacement  of  the  tip  into  the  sample,  and 
S*=  S.  The  contact  depth,  hc ,  is  related  to  the  deforma¬ 
tion  behavior  of  the  material  and  the  shape  of  the 
indenter,  as  illustrated  in  Fig.  7,  and  is  given  by 
K  -  ^max-  K  >  where  hs  is  defined  as  the  elastic  displace¬ 
ment,  sometimes  referred  to  as  sink-in,  of  the  surface  at 
the  contact  perimeter.  For  each  of  three  specific  tip 
shapes  (flat-ended  punch,  paraboloid  of  revolution,  and 
cone),  hs  =  ePmax/S  where  £  is  a  function  of  the  partic¬ 
ular  tip  geometry,  as  summarized  in  Table  1.  Thus,  hc  is 
given  by 


Despite  the  use  of  feedback  in  some  cases,  the 
dynamic  model  of  Fig.  5  is  normally  treated  as  a 
simple  damped  harmonic  oscillator  under  conditions 
of  an  applied  harmonic  force,  P  =  P0-  exp  ( icot ),  with 
a  resulting  harmonic  displacement,  h  =  h0-  exp  ( icot -  0), 
where  mis  the  applied  oscillation  frequency  and  0is  the 
phase  angle  by  which  the  displacement  lags  the  force. 
The  equations  of  motion  developed  for  this  model  (with 
sample  contact)  can  be  solved  to  yield  an  equation  for 
the  contact  stiffness,  S ,  that  is  a  function  of  previously 
calibrated  instrument  parameters  (i.e.,  Kf,  Ks,  and  m), 
co,  0,  the  harmonic  force  amplitude,  P0,  and  the  har¬ 
monic  displacement  amplitude,  h0: 


S  = 


—cos  0  -  (Ks  -  mco2 ) 

K 


Kf 


(7) 


Also,  the  damping  factor,  C,  attributed  to  the  tip- 
sample  contact  is  given  by: 


Cco  =  —  sin  0- Dm. 

K 


(8) 


For  the  purposes  of  the  Oliver-Pharr  analysis,  h  -  hmax 
and  P  =  Pmax  in  Eq.  (6).  Also,  hc  <  hmax  such  that  this 
equation  is  not  valid  when  pile-up  occurs,  i.e.,  when 
material  is  forced  up  along  the  sides  of  the  indentation 
tip.  As  indicated  in  Table  1,  the  choice  of  £  should  be 
related  to  the  value  of  m  determined  from  the  curve 
fitting  [40].  However,  a  value  of  0.75,  corresponding 
to  m  =  1.5,  is  used  almost  exclusively  for  spherical, 
conical,  and  pyramidal  indenters,  regardless  of  the 
value  of  m.  Once  hc  has  been  determined,  the  tip  shape 
function,  A(hc)  (see  Sec.  2.4)  is  used  to  calculate  the 
contact  area,  such  that  the  sample  modulus,  E ,  can  be 
determined  from  Eq.  (3)  and  Eq.  (4),  given  a  reasonable 
estimate  of  v. 


Superposing  oscillation  during  quasi- static  loading  thus 
allows  S  to  be  estimated  as  a  function  of  depth  through¬ 
out  the  loading  segment.  Equation  (6)  can  then  be  used 
to  monitor  hc  continuously,  knowledge  of  the  tip  shape 
yields  an  estimate  of  A  continuously,  and  thus  E  can  be 
measured  continuously  using  Eq.  (3)  and  Eq.  (4). 
Because  these  relationships  can  be  built  into  the  soft¬ 
ware,  the  use  of  dynamic  oscillation  superposed  over 
the  loading  history  simplifies  the  analysis  procedures 
greatly.  Additionally,  the  amount  of  data  associated 
with  one  such  experiment  is  equivalent  to  many 
experiments  in  which  only  quasi- static  loading  is  used, 
allowing  more  robust  statistical  analysis  to  be 
performed. 


Fig.  7.  Illustration  of  the  indentation  geometry  at  maximum  force  for  an  ideal 
conical  indenter. 
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2.4  Calibration  Issues 

Very  little  discussion  of  calibration  issues  related  to 
instrumented  indentation  systems  can  be  found  in  the 
open  literature.  As  with  most  measurement  systems, 
calibration  is  essential  for  limiting  uncertainties  and 
achieving  reproducible  and  repeatable  measurements. 
The  most  fundamental  measures  made  by  these  systems 
are  of  force  and  displacement.  Many  instruments  are 
capable  of  operating  at  forces  less  than  1  mN,  and  in 
some  literature,  force  resolutions  of  1  pN  or  less  are 
claimed.  Currently,  however,  force  is  typically  calibrat¬ 
ed  by  hanging  standard  weights  on  the  force  measure¬ 
ment  device,  and  deadweight  force  standards  are  only 
available  for  calibrating  forces  down  to  approximately 
10  pN.  Current  NIST  research  aimed  at  improving 
available  force  standards  at  the  micro-Newton  and 
nano-Newton  levels  is  discussed  in  Sec.  3.1. 

Calibration  of  displacement  can  be  done  using  sever¬ 
al  methods,  such  as  by  using  separately  calibrated 
transducers  or  by  using  interferometric  methods. 
Again,  however,  many  instruments  are  capable  of  oper¬ 
ating  with  displacements  and  displacement  resolutions 
significantly  smaller  than  the  resolution  of  the  calibra¬ 
tion  methods.  Further,  force  and  displacement  measure¬ 
ments  are  typically  coupled  in  instrumented  indentation 
systems  through  the  support  springs.  Thus  for  force- 
driven  systems,  for  example,  calibrating  displacement 
is  equivalent  to  calibrating  the  spring  stiffness  or  spring 
constant,  which  links  the  raw  displacement  of  the 
system  to  the  raw  force.  Assuming  the  spring  response 
is  linear,  displacements  can  then  be  determined  with 
uncertainties  related  to  uncertainties  in  the  force 
calibration  and  the  spring  constant. 

The  calibration  of  load-frame  compliance,  Clf,  is 
difficult  and  can  have  a  significant  amount  of  uncer¬ 
tainty,  and  even  qualitative  assessment  of  indentation 
behavior  depends  critically  on  the  accuracy  of  this  cal¬ 
ibration  step,  especially  from  laboratory  to  laboratory 
and  instrument  to  instrument.  Prior  to  calibration  (i.e., 
Clf  is  unknown),  the  measured  displacement,  /ztotal,  is  a 
combination  of  displacement  of  the  load  frame,  /%,  and 
displacement  of  the  sample,  hsp .  Treating  the  system  as 
two  springs  (the  load  frame  and  the  sample)  in  series 
under  a  given  force,  P , 

Ki*=h+K-  (9) 

Dividing  both  sides  by  P , 

C“al  =  C]i  +  C*  =  ^  +  2E;  7T  (10) 


The  total  compliance  is  related  to  total  stiffness,  S*9  by 
Ctotai  =  1  /S*  and  the  sample  compliance  is  related  to  the 
contact  stiffness,  S ,  by  Csp  =  1/S.  A  number  of  possible 
methods  exist  for  determining  Clf  using  reference 
samples  that  are  homogeneous  and  isotropic  and  for 
which  both  E  and  v  are  known.  Typically,  a  series  of 
indentation  measurements  are  made  on  a  single  refer¬ 
ence  sample  or  multiple  reference  samples.  Oliver  and 
Pharr  [28]  suggested  using  an  iterative  technique  to 
calibrate  both  the  load-frame  compliance  and  the  tip 
shape  with  one  set  of  data  from  a  single  reference 
sample,  as  both  Clf  and  A  are,  in  general,  unknowns  in 
Eq.  (10).  While  this  method  has  the  advantage  of  not 
requiring  an  independent  measurement  of  the  area  of 
each  indent,  its  use  has  been  limited,  perhaps  because  it 
is  mathematically  and  time  intensive. 

For  the  load-frame  compliance  calibration,  the 
choice  of  reference  sample(s)  should  be  made  with  an 
objective  of  maximizing  contact  stiffness  (minimizing 
Csp)  so  that  Ctotal  is  dominated  by  Clf.  Thus,  relatively 
large  indentation  forces  and  depths  are  normally 
applied  to  a  reference  material  that  either  has  a  high 
modulus  or  exhibits  significant  plastic  deformation, 
i.e.,  a  sample  that  has  a  high  ratio  of  E/H.  Use  of  a  high¬ 
ly  plastic  material  such  as  aluminum,  for  which  the  pro¬ 
jected  area  should  be  similar  to  the  contact  area,  A,  at 
maximum  force,  can  be  combined  with  high  resolution 
imaging  techniques,  such  as  electron  microscopy  or 
atomic  force  microscopy  (AFM)  to  determine  Clf  from 
a  plot  of  Ctotal  as  a  function  of  1  /J~A  .  Alternatively,  an 
additional  assumption  that  H  is  constant  with  depth 
can  allow  Clf  to  be  determined  from  a  plot  of  Ctotal  as  a 
function  of  1  lJPmax  •  In  this  case,  aluminum  is  often 
replaced  by  fused  silica,  because  oxide  formation  on 
aluminum  can  create  variations  in  both  E  and  H  with 
penetration  depth.  For  all  of  these  methods,  the  calcu¬ 
lation  of  Clf  based  on  Eq.  (10)  requires  a  large  extrapo¬ 
lation  of  the  experimental  data.  Futher,  the  v- variable 
has  significant  uncertainty,  which  is  a  violation  of  the 
assumptions  of  least  squares  regression,  leading  to  the 
creation  of  bias  in  the  least  squares  estimate  of  Clf  such 
that  the  estimated  value  is  lower  than  the  actual  value. 
These  issues,  along  with  other  sources  of  uncertainty, 
can  result  in  a  large  amount  of  uncertainty  associated 
with  the  load-frame  compliance,  which  will  then  affect 
all  subsequent  calculations  [20,  41].  A  review  of  these 
uncertainties  along  with  a  newly  proposed  method  of 
load-frame  compliance  calibration  that  is  independent 
of  reference  materials  can  be  found  in  Ref.  [42]. 

For  the  tip  shape  calibration,  the  series  of  indents 
applied  to  a  reference  material  typically  covers  a  larger 
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range  of  maximum  force  and  maximum  penetration 
depth  compared  to  the  load-frame  compliance  calibra¬ 
tion  procedure.  The  objective  of  tip  shape  calibration  is 
to  estimate  the  cross-section  area  of  the  indenter  tip  as 
a  function  of  distance  from  the  apex.  In  Fig.  7,  the 
indentation  geometry  for  a  conical  indenter  is  illustrat¬ 
ed  in  two  dimensions.  At  a  given  force,  P,  the  contact 
area,  A,  which  is  related  to  the  contact  radius,  a ,  is  the 
cross-section  area  of  the  indenter  tip  at  a  distance,  hc 
(the  contact  depth),  from  the  tip  apex.  From  measure¬ 
ments  of  /zmax,  Pmax,  and  S ,  Eq.  (3)  and  Eq.  (6)  can  be 
used  to  calculate  A  and  hc ,  respectively,  for  each  inden¬ 
tation.  A  tip  shape  function,  A(/zc),  is  determined,  given 
a  sufficient  number  of  measurements  over  a  range  of  hc 
values,  by  fitting  the  A  vs  hc  data,  typically  using  a 
multi- term  polynomial  fit  of  the  form: 

A(hc)=BX  +  Blhc  +  B2hl,2+....  (11) 

In  this  equation,  B0 ,  Bn  are  constant  coefficients 
determined  by  the  curve  fit.  For  example,  for  a 
Berkovich  indentation  tip,  Oliver  and  Pharr  [28]  sug¬ 
gested  using  up  to  nine  terms  (n  =  8)  with  B0  =  24.5, 
where  the  area  function  of  a  perfect  Berkovich  tip  is 
A{hc )  =  24.5  hc 2.  The  additional  terms  attempt  to 
account  for  deviations  from  ideal  geometry,  such  as 
blunting  of  the  tip. 

The  use  of  dynamic  oscillation  can  significantly 
enhance  the  calibrations  of  load-frame  compliance  and 
tip  shape,  particularly  with  regard  to  improving  the 
statistical  sampling  and  reducing  the  amount  of  time 
required.  For  example,  multiple  deep  indentations  on  a 
reference  sample  using  oscillation  superposed  over  the 
loading  segment  of  each  test  yields  multiple  sets  of 
modulus  values  as  a  function  of  penetration  depth  that 
can  then  be  averaged  and  used  either  to  check  load- 
frame  compliance  or  tip  shape  calibration.  In  fact,  both 
can  at  least  be  checked  with  the  same  data  if  the  refer¬ 
ence  sample  (e.g.,  fused  silica  glass)  has  both  a  modu¬ 
lus  and  a  hardness  that  is  independent  of  depth.  Thus, 
the  E  vs  h  data  can  be  used  to  check  tip  shape  calibra¬ 
tion,  and  for  a  proper  load-frame  compliance  cali¬ 
bration,  the  ratio  of  force  to  contact  stiffness  squared 
(P/52),  which  is  proportional  to  H/Er ,  should  be  inde¬ 
pendent  of  depth  regardless  of  the  tip  shape  function. 
Additionally,  the  improved  capabilities  for  detecting 
the  surface  using  dynamic  oscillation  reduce  the  uncer¬ 
tainties  in  the  calibrations  related  to  the  choice  of  the 
initial  point  of  contact. 

Dynamic  calibrations  of  the  system  are  typically 
made  with  respect  to  the  dynamic  model  shown  in 


Fig.  5  [37].  System  calibrations  are  then  performed  by 
measuring  the  dynamic  response  with  no  sample 
involved.  The  load-frame  stiffness,  Kf,  is  normally 
assumed  to  be  infinite  (i.e.,  load- frame  compliance  is 
zero).  By  monitoring  amplitude  and  phase  shift,  the 
equations  derived  for  the  model  can  be  used  to  deter¬ 
mine  the  resonance  frequency  of  the  system,  the  system 
damping  coefficient,  D ,  the  mass,  and  the  spring 
constant,  Ks .  As  discussed  previously,  the  spring  con¬ 
stant,  which  is  typically  independent  of  frequency  over 
a  wide  frequency  range,  can  be  determined  from  raw 
force  as  a  function  of  raw  displacement,  and  the  system 
mass  can  be  determined  from  the  displacement  at  zero 
force.  Also,  D  is  often  assumed  to  be  independent  of 
frequency,  which  is  not  necessarily  true  [43]. 

3.  Research  Efforts  at  NIST 

3.1  Improving  Force  Calibration 

A  desire  for  accurate,  traceable,  small  force  measure¬ 
ment  is  emerging  within  the  International  Organization 
for  Standardization  (ISO)  task  groups  and  ASTM 
International  technical  committees  that  work  on  instru¬ 
mented  indentation  standards  [44] .  However,  no  meth¬ 
ods  for  establishing  force  measurement  traceability  at 
these  levels  are  currently  available.  A  research  effort 
has  recently  been  developed  at  NIST  (Microforce 
Realization  Competence)  with  the  purpose  of  creating  a 
facility  and  instruments  capable  of  providing  a  viable 
primary  force  standard  below  1 0“5  N,  and  with  the  goal 
of  realizing  force  in  this  range  at  a  relative  uncertainty 
of  parts  in  104.  This  new  project  complements  a  body  of 
existing  work  at  NIST  to  develop  standards  and  meth¬ 
ods  for  the  instrumented  indentation  community  that 
together  provide  a  metrological  basis  for  manufacturers 
seeking  traceable  characterization  of,  for  example,  thin 
film  mechanical  properties. 

The  most  common  approach  to  force  realization  is  a 
calibrated  mass  in  a  known  gravitational  field  or  dead¬ 
weight  force,  which  is  universally  accepted  as  the 
primary  standard  of  force.  The  smallest  calibrated 
mass  available  from  NIST  is  1  mg  (approximately  a 
10  pN  deadweight  force)  having  a  relative  uncertainty 
of  about  10^.  In  principle,  smaller  masses  could 
be  calibrated,  but  they  would  be  difficult  to  handle. 
Also,  the  relative  uncertainty  tends  to  increase  in 
verse  proportion  to  the  decrease  in  mass  [45],  potential¬ 
ly  resulting  in  uncertainties  that  are  of  similar  magni¬ 
tude  to  deadweight  forces  in  the  range  of  nano- 
Newtons. 
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Alternatively,  forces  in  this  range  can  be  realized 
using  the  electrical  units  defined  in  the  International 
System  of  Units  (SI)  and  linked  to  the  Josephson 
and  quantized  Hall  effects  in  combination  with  the 
SI  unit  of  length.  This  realization  can  be  done  using 
electromagnetic  forces  (e.g.,  the  NIST  Watt  Balance 
Experiment  [46])  or  using  electrostatic  forces  [47].  In 
this  research,  the  latter  was  chosen  because  the  required 
metrology  is  somewhat  simpler  to  execute,  and  the 
forces  generated,  although  generally  less  than  those 
feasible  electromagnetically,  are  appropriate  for  the 
force  range  of  interest.  Also,  electrostatic  force  genera¬ 
tion  is  common  in  micro-electromechanical  systems 
(MEMS),  and  the  ability  to  calibrate  such  forces 
from  electrical  and  length  measurements  could  prove 
beneficial. 

The  mechanical  work  required  to  change  either  the 
overlap  or  the  separation  of  two  electrodes  in  a  one¬ 
dimensional  capacitor  while  maintaining  constant 
voltage  is 

dW  =  F  dz  =  I V2dC .  (12) 

In  this  equation,  d W  is  the  change  in  energy  (mechani¬ 
cal  work),  F  is  the  force,  d z  is  the  change  in  the  overlap 


or  separation  of  the  electrodes,  V  is  the  electric  poten¬ 
tial  across  the  capacitor,  and  dC  is  the  change  in  capac¬ 
itance.  Thus,  force  can  be  realized  from  electrical  units 
by  measuring  V  and  the  capacitance  gradient,  dC/dz  or: 


This  idealized,  one-dimensional  approximation  does 
not  account  for  multi-dimensionality,  external  fields  or 
stray  electrical  charges  likely  to  be  present  in  the 
actual  physical  system.  The  goal  then  is  to  develop  a 
system  that  reproduces  this  idealization  as  closely  as 
possible  by  using  effective  constraints  on  the  geometry, 
shielding,  and  suspension  of  the  resulting  electrodes. 
Additionally,  validation  of  this  electrostatic  force 
realization  through  comparison  to  deadweight  forces 
will  be  advantageous,  at  least  in  the  higher  force  range 
where  the  uncertainty  that  can  be  achieved  mechanical¬ 
ly  is  still  competitive.  In  consideration  of  these  factors, 
a  force  generator  was  designed  to  operate  along  the 
vertical  axis  as  part  of  an  electromechanical  null 
balance  shown  in  Fig.  8.  Recent  results  with  this  instru¬ 
ment  [48,  49],  referred  to  as  the  NIST  Electrostatic 
Force  Balance  or  EFB,  demonstrated  a  relative 


Fig.  8.  The  prototype  electrostatic  force  balance.  Inner  cylindrical  electrode  of  15  mm  dia.  is 
suspended  from  a  compound  parallelogram  leaf  spring  made  of  50  pm  thick  CuBe  producing  a 
single  axis  spring  of  stiffness  13.4  N/m.  Deflections  are  measured  using  a  double-pass  Michelson 
interferometer  and  nulled  using  a  feedback  servo  to  apply  voltage  to  the  outer  cylinder.  Electrode 
gap  is  nominally  0.5  mm  and  overlap  is  nominally  5  mm. 
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standard  uncertainty  of  about  10"4  in  the  comparison  of 
gravitational  and  electrostatic  forces  ranging  between 
10  pN  and  100  pN.  This  result  indicates  that  the  elec¬ 
trostatic  force  can  be  constrained  and  measured  in  a 
fashion  traceable  to  the  SI  and  with  accuracy  sufficient 
to  warrant  consideration  as  a  primary  standard  of  force 
in  this  regime. 

3.2  Improving  Tip  Shape  Calibration 

A  number  of  recent  efforts  have  been  made  to 
improve  tip  shape  calibration  for  instrumented  indenta¬ 
tion  [41,  50-55].  These  efforts  have  included  material- 
independent  methods  of  tip  shape  calibration  using 
AFM  [41,  50-52]  and  alternative  procedures  using 
indentation  of  reference  materials  [53-55].  Recent 
research  at  NIST  has  focused  on  methods  in  which  the 
indenter  tip  is  scanned  with  an  AFM  probe  to  yield 
direct  information  regarding  the  three-dimensional  tip 
shape  [53,  56].  Because  an  AFM  image  is  a  combina¬ 
tion  of  the  AFM  probe  geometry  and  the  geometry  of 
the  sample  surface,  AFM  imaging  of  indentation  tips 
was  combined  with  AFM  imaging  of  a  tip  characteriz- 
er  surface,  which  allows  for  an  estimation  of  the  three- 
dimensional  shape  of  the  AFM  probe  using  the  method 
of  blind  reconstruction  [57,  58].  This  estimation  of  the 
AFM  probe  geometry  can  then  be  eroded  from  the 
image  of  the  indentation  tip  to  remove  dilation  and 
other  artifacts  created  by  the  AFM  probe,  as  described 
in  detail  in  Ref.  [50]. 


In  principle,  because  the  AFM  probe  acts  to  produce 
image  geometry  that  is  dilated  from  the  true  surface 
geometry,  the  area  function  (cross-section  area  as  a 
function  of  distance  from  the  apex)  of  the  indentation 
tip  determined  from  an  AFM  image  is  an  outer  bound 
on  the  true  area  function.  Also,  because  the  AFM  probe 
geometry  estimated  using  blind  reconstruction  is  an 
outer  bound  on  the  true  probe  geometry,  eroding  this 
estimated  geometry  from  an  AFM  image  of  the  inden¬ 
tation  tip  produces  a  lower  bound  on  the  tip  area  func¬ 
tion.  In  reality,  however,  other  artifacts  or  uncertainties 
associated  with  the  AFM  image  of  the  indentation  tip 
will  affect  the  results,  particularly  for  open-loop  AFM 
systems.  This  effect  is  shown  in  Fig.  9  in  which  area 
functions  are  compared  for  consecutive  images  of  the 
same  size  and  of  different  sizes.  For  close-loop  AFM 
systems,  particularly  those  with  calibrated  vertical 
motion  such  as  the  NIST  Calibrated- AFM  [59],  the 
images  produced  of  the  tip  shapes  can  have  much  less 
uncertainty. 

Comparisons  of  tip  shape  data  generated  from  AFM 
imaging  to  that  determined  from  indentation  of  fused 
silica  (for  example,  see  Fig.  9)  revealed  a  potential 
problem  with  determining  tip  shape  area  functions 
using  indentation  of  reference  samples.  Differences 
between  AFM-generated  data  and  indentation  data 
were  most  significant  at  small  depths  and  increased  as 
a  function  of  the  indenter  tip  radius.  Values  of  cross- 
section  area  from  indentation  data  were  always  less 
than  that  from  AFM  data.  This  result  appears  to  at  least 
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Fig.  9.  Plot  of  the  estimated  cross-section  area  of  a  Berkovich  indentation  tip  as  a  function  of  the  distance  from  the  apex.  Data  are 
shown  for  two  3  pm  x  3  pm  images,  one  16  pm  x  16  pm  image,  each  of  the  three  AFM  images  after  eroding  away  the  estimated  AFM 
probe  geometry,  and  indentation  of  fused  silica.  Only  the  first  300  nm  of  data  are  shown  to  emphasize  differences. 
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partially  explain  results  in  which  indentation  modulus 
(E  oc  SI  Va)  is  significantly  higher  than  modulus  values 
measured  with  other  techniques,  particularly  for  poly¬ 
mers  where  S  is  relatively  small  for  a  given  contact 
depth  such  that  errors  in  A  have  a  larger  effect  and  low 
values  of  A  will  result  in  artificially  high  modulus 
values.  Also,  this  results  could  also  explain  the  many 
reports  of  indentation  hardness  ( H  =  P/A)  values  being 
larger  at  small  depths  compared  to  large  depths,  often 
referred  to  as  the  indentation  size  effect.  Thus,  the 
uncertainties  associated  with  tip  shape  calibration  using 
indentation  of  reference  samples  are  expected  to  be 
significantly  larger  than  those  associated  with  AFM- 
generated  tip  shape  information  in  many  cases.  In 
fact,  this  independent  assessment  of  the  tip  shape  via 
AFM  can  then  be  used  to  help  calibrate  load-frame 
compliance  using  a  known  tip  area  function. 

3.3  Applications  of  Instrumented  Indentation  to 
Polymeric  Materials 

Instrumented  indentation  is  increasingly  being  used 
to  probe  the  mechanical  response  of  polymeric  and 
biological  materials.  These  types  of  materials  behave  in 
a  viscoelastic  fashion,  i.e.,  they  display  mechanical 
properties  intermediate  between  those  of  an  elastic 
solid  and  a  viscous  fluid,  as  shown  in  Fig.  2  and  Fig.  3 
by  the  changes  in  displacement  under  constant  force 
conditions  and  in  Fig.  4  by  the  changes  in  force  under 
constant  displacement  conditions.  The  mechanical 
behavior  is  thus  dependent  on  the  test  conditions, 
including  amount  of  strain,  strain  rate,  and  temperature. 
Often  in  instrumented  indentation,  however,  properties 
are  measured  using  loading  histories  developed  for 
elastic  and  elasto-plastic  materials,  the  properties  of 
which  are  not  particularly  time  dependent,  and  the 
analysis  of  the  indentation  response  is  typically  based 
on  elasticity  theory.  In  studies  in  which  attempts  have 
been  made  to  characterize  viscoelastic  behavior,  limit¬ 
ing  and  sometimes  invalid  assumptions  have  been 
made,  and  linear  viscoelasticity  has  been  applied 
despite  the  intense  strains  local  to  the  indenter  tip 
that  would  appear  to  violate  the  linear  viscoelasticity 
premise  of  infinitesimal  strains  [60] . 

Another  issue  that  can  add  significant  uncertainty  to 
indentation  measurements  of  polymers  is  the  ability  to 
detect  the  surface.  Polymers  and  other  organic  materi¬ 
als  are  typically  much  more  compliant  compared  to 
the  metallic  and  ceramic  types  of  materials  to  which 
instrumented  indentation  mainly  has  been  applied,  with 
modulus  values  ranging  from  a  few  GPa  for  common 


glassy  polymers  to  a  few  MPa  or  lower  for  rubbery 
polymers  and  many  biological  materials.  As  mentioned 
previously,  the  use  of  small  dynamic  oscillations  has 
improved  sensitivity  to  surface  contact.  However,  as 
the  compliance  of  the  material  increases,  depending  on 
the  particular  indentation  system,  the  changes  used  to 
define  surface  contact  become  of  similar  magnitude  to 
the  noise  level.  In  some  cases,  manual  selection  of  the 
contact  point  in  the  raw  test  data  after  the  test  is  com¬ 
pleted  can  be  used  to  correct  any  automated  selection 
by  the  system  or  to  check  the  sensitivity  of  the  calcula¬ 
tions  to  this  choice.  However,  in  other  cases,  the  system 
selection  of  the  contact  point  affects  the  start  of  the  test 
and  thus  the  start  of  any  feedback  that  might  be  used 
to  control  the  test  flow.  In  such  cases,  the  approach 
velocity  of  the  probe  can  be  an  important  factor,  as  the 
initial  part  of  the  indentation  data  prior  to  feedback 
control  will  be  based  on  this  approach  velocity.  For 
example,  the  indentation  strain  rate,  ,  as  estimated 
by  the  ratio  PIP  ,  and  modulus  are  plotted  in  Fig.  10 
as  functions  of  penetration  depth  for  a  polystyrene 
material.  The  feedback  used  to  keep  PIP  constant 
does  not  take  effect  until  the  probe  is  over  100  nm  into 
the  material.  Although  other  factors  contribute  to 
significant  uncertainty  at  small  depths,  the  large 
changes  in  indentation  strain  rate  could  also  have 
affected  the  resulting  modulus  values. 

In  recent  research  at  NIST,  analyses  by  Ting  [61]  that 
are  based  on  contact  between  a  rigid  indenter  and  a 
linear  viscoelastic  material  were  revisited  and  used  to 
determine  under  what  conditions,  if  any,  instrumented 
indentation  can  be  used  to  measure  linear  viscoelastic 
behavior  for  a  number  of  different  polymers  [58].  For 
example,  the  creep  compliance,  J{t\  of  a  linear  vis¬ 
coelastic  material  subject  to  a  constant  indentation  or 
creep  force,  P0,  using  a  conical  tip  of  semi-apical  angle, 
6 ,  is  proportional  to  the  change  in  contact  area  with 
time,  A(t)  [61]: 


J(f)  oc 


A(t) 
P0  tan  6 


(14) 


Note  that  in  this  equation,  (1/tan  6)  is  related  to  a  nom¬ 
inal  indentation  strain  and  P/A  is  related  to  a  nominal 
indentation  stress,  such  that  J{t)  is  proportional  to  an 
indentation  strain  over  an  indentation  stress.  Similar 
equations  can  be  determined  for  other  tip  geometries, 
and  equations  can  be  derived  relating  stress  relaxation 
modulus  to  the  change  in  force  with  time  during  a 
constant  displacement  indentation  test  for  various  tip 
geometries. 
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Fig.  10.  Plots  of  indentation  strain  rate,  e. ,  estimated  by  the  ratio  PIP  (a)  and  indentation  modulus 
(b)  as  a  function  of  depth,  h,  for  a  constant  p/p  test  using  a  Berkovich  indentation  tip  to  pene¬ 
trate  a  polystyrene  sample.  Date  for  one  test  each  at  two  different  rates  are  shown.  Differences  in  the 
modulus  data  are  estimated  to  be  within  the  measurement  uncertainty. 


In  these  studies,  constant  force  indentation  creep 
tests  and  constant  penetration  depth  stress  relaxation 
tests  were  used,  and  the  results  were  compared  to 
traditional  solid  rheological  measurements  [56].  An 
example  of  indentation  creep  compliance  measure¬ 
ments  for  an  epoxy  material  using  a  rounded  conical 
tip  (manufacturer-determined  tip  radius  of  10  pm)  is 
shown  in  Fig.  11  [proportionality  constant  of  1.0 
assumed  in  Eq.  (14)].  The  dependence  of  creep  compli¬ 
ance  on  the  creep  force,  as  shown  by  the  displacement 
between  the  curves,  is  an  indication  of  nonlinear  behav¬ 


ior,  and  in  fact,  the  indentation  creep  behavior  of  a 
number  of  polymeric  materials  was  dominated  by 
nonlinear  viscoelastic  behavior.  Additionally,  probe  tip 
size  and  shape  were  altered  to  produce  different  nomi¬ 
nal  indentation  strains,  e.  ,  and  the  measured  respons¬ 
es  appeared  to  be  correlated  with  £  .  Analyses  and 
protocols  are  currently  being  explored  for  relating 
instrumented  indentation  data  to  viscoelasticity  and 
ultimately  to  stress-strain  behavior  of  polymers  using 
various  indentation  tips. 
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Fig.  11.  Log-log  plot  of  creep  compliance,  J(t),  as  a  function  of  time,  t,  for  an  indentation  creep  experiment  on  epoxy  using  a  rounded 
conical  tip  (manufacturer-determined  tip  radius  of  10  pm).  Error  bars  shown  represent  an  estimated  standard  deviation  (k  =  1). 


In  addition  to  the  quasi- static  studies,  measurements 
of  viscoelastic  behavior  using  dynamic  indentation 
techniques  are  also  being  explored.  From  the  dynamic 
model  of  the  indentation  system  (see  Fig.  5),  equations 
have  been  derived  for  determining  the  storage  modulus, 
E\  and  loss  modulus,  E'\  of  a  viscoelastic  material  [18, 
19,  37]: 

/ 

E^S,  (15) 


E'  = 


(16) 


Thus,  E'  is  assumed  to  be  directly  related  to  the  storage 
portion,  S ,  of  the  mechanical  impedance  of  the  tip- 
sample  contact  in  the  dynamic  model  (see  Fig.  5),  and 
E"  is  assumed  to  be  directly  related  to  the  loss  portion, 
Cco ,  of  the  mechanical  impedance  of  the  tip- sample 
contact.  However,  this  assumption  was  found  to  have 
significant  limitations  with  regard  to  lossy  polymers 
[43].  At  NIST,  new  analysis  methods  are  currently 
being  developed  for  analyzing  the  dynamic  mechanical 
response  to  indentation  of  viscoelastic  materials,  in 
particular  to  determine  under  what  conditions  Eq.  (15) 
and  Eq.  (16)  hold.  A  number  of  different  polymers  are 
being  studied,  the  results  of  which  are  being  compared 
to  traditional  dynamic  mechanical  measurements. 


4.  Summary 

As  instrumented  indentation  techniques  gain  wider 
use  and  application,  standardization  efforts  increase  in 
importance.  NIST  personnel  are  involved  in  the  ASTM 
Task  Group  E28.06.ll  developing  ASTM  standard 
practices  and  standard  test  methods  for  instrumented 
indentation  testing  and  have  participated  in  internation¬ 
al  round  robin  testing.  Deficiencies  in  current  practices 
include  the  need  to  use  reference  materials  for  calibra¬ 
tions  of  load-frame  compliance  and  tip  shape,  the  lack 
of  traceable  force  calibration  below  10  pN,  and  the  lack 
of  uncertainty  budget  analysis  related  to  measurement 
techniques  and  analyses.  Current  research  at  NIST  is 
focused  on  independent  methods  for  tip  shape  calibra¬ 
tion,  traceable  calibration  methods  for  micro-Newton 
and  nano-Newton  level  forces,  and  applications  to 
viscoelastic  materials  such  as  polymeric  and  biological 
materials. 
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